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A^TFy\CT 

Computer Aided Tonography (CAT) is vddely used in 
the medical imaging area for tte detection of cancerous 
tissues. It has also proved to be a good technique for 
measuring point-density (void-fraction) in two-phase flow 
systons. The present da y GAT-scanners employ algorithms 
which are essentially a discrete implementation of the Radon 
Inversion Formulae, These methods of reconstructing density 
distributions are fast and easy to implement. 

The present work is an attempt to simplify (and conse- 
quently to speed-up) the existing algorithms for the cases of 
radially symmetric density distribution. The method uses 
Radon Inversion Formulae and involves Bessel’s Function of 
Fourier frequency. 

The algorithm was tested on some known distributions 
and the reconstruction is excellent. Data for bubbly-flow was 
also processed to obtain a point density distribution. The 
range of investigation wbs 0,6^ <P> < 0*9 gm/cm , wtiere 
" < > ” denotes average cross-sectional value. The values 
of <Pqj> are within + 0,03 g/cm of the true values, <p>, 
for the data sets coreidered. The results a3re comparable with 
another simple method involving radial polyndmials. It is 
observed that the simplified Radon Inversion Formulae for 
the radially syrame trie case gives good results and can be 
implemented with ease, , 



1 


CHAPTER 1 
INTRDDUCTION 

- . nk Ik % ,m Hi.** . * I — 

The problem of image reconstruction from projections 
has repeatedly arisen over the last four decades in a large 
number of scientific, medical and technical fields. Of all 
the applications, probably the greatest effect on the world 
at large has been in the area of diagnostic medicine^ 
Gompute 2 >-Aided-Tomography (GAT) has revolutionized tJ« 
techniques in radiology. 

Images of cross-sections of the human body are 
produced from the data obtained by measuring the attenuation 
of radiation along a large number of lines through the 
cross-section. The CAT machine incorporates scanning of 
patient with radiation using appropriate tomographic algorithm 
to reconstruct the density distribution of tissues [1] , 

A major problem confronting an engineer in the 
Nuclear reactor safety field is the measurement of void 
fraction (or density) during various accident/ transient 
situations. The developHnaent of an accurate and ireliable 
two-i^ase flow measurement technique for various flow systenB 
and ccwiponents in Light Water Reactor (L\®.) is of great 
importance and significance. The concept of measuring ctensity 
distribution was first investicpited by Schlosser et ai [2] 
for a two-j^ase air-water system. The density measurement 



results have been summarized by Kulacki et al [ 3 

In many applications where images are required, 
the n^asurements can be done by probing tlte object by 
invisible radiations and interpreting the measured data. 
Sometimes immediate interpretation of measured data is not 
possible directly, but is related to the image in a known 
way. Generally, all image reconstruction procedures or 
algorithms aim at processing the data to form the image and 
thereby to facilitate the interpretation of measurement. 

The present work mainly looks at such a new field, 
which will alleviate the problams (related to radially 
symmetric distribution) in the nuclear engineering area (two- 
phase flow in reactor systems), chemical industry, food 
processing etc. 

The most important feature of tomographic recons- 
truction method is its non-invasive aspect of polpt-density 
measurement. Due to this fact, the system retains its 
natural/ original configuration, without being even slightly 
affected by the measuring device. 

The present work is aimed at developing a simplified 
algorithm based on the Radon Inversion Formulae, involving 
Bessel Functions, which can perform tomographic reconstruction 
of the density distribution for radially symmetric two-ihase 
flow®. The data is collected in the parallel-beam mode. 



The algorithm written in FCSITRAN^ has been initially tested 
against some known functions, simulated radially symmetric 
distributions, and annular flows. The algorithm has also 
been successfully implemented to reconstruct the density 
map for air-water bubbly flow data [3,4] , The results 
are comparable with the results obtained by the other 
known general methods [3,4,53 • The results are also 
comparable with another alternate scheme involving radial 
polynomials [ 6,7] , 
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CHAPTER 2 

— A •mu* JT -■fc.tA x* Jt 

THEORETICAL FORMULATION 

- w -m J0. a x:i%. ■« Mr. jrtx u 

This chapter includes a brief discussion of: 

(a) methods of data collection in a tomographic 
system depending upon geometry and, 

(b) the relation of this data to the distribution. 

Two standard geometries of data collection, namely, 
fan beam and parallel beam geometries are discussed. In 
addition, types of algorithms, two phase flow and use of 
CAT in identifying flow regimes are briefly explained, 

2,1 Data Collection 

In a tomographic system, radiation is passed through 
the object (fo be -lleconstructed) and the attenuated signal 
is detected. This signal is related to the function to be 
reconstructed. The single beam radiation attenuation 
phenomenon is represented by: 

N = N^ exp t- ^ |i(r,0)dsl (1) 

where, 

N = detector reading (counts/s) 

Nq ss source strength (counts/s) 

S = path of radiation (ray) 
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C = chord along which S is integrated 
}i = absorption coefficient 
r,j& = cylindrical co-ordinates. 

Rewriting eqn ,( 1 ), we have 

p = ^ |i(r,0)ds (2) 

where 

p = In (Nq/N) (3) 


2.2 


Parallel Beam Geometry Data Collection 


A line in the plane is specified by two parameters; 
its (signed) distance 1 from the origin and its angle © 
with jrespect to the Y-^xis (Fig.l). The function of tvro 
variables, whose value for any (1,0) is defined as line- 
integral of p, along the line specified by 1 and 0, is denoted 
by p=(l,©). For any 1 and ©, (1,0) and (-1, 8+7t) represents 
the same line in the plane, so that 


p(l,0) = p(-l, ©+15 ) 


(4) 


By projection theorem [8] we know that 
^(R Cos©, R Sin©) = p(R,0) 

where 

^R,©) » ? p( l,©) e~^^ dl 

-D 


(5) 

( 6 ) 
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PARALLEL BEAM GEOMETRY 



LEGEND: S - Source 

D - Detector 

Points within object arc described by a fixed 
(x,y)or (r co-ordinate system- Rays 

(dashed-line) are specified by their angle with 
respect toy axis,8 and their distance from origin 1. 


Fig.1 Data collection geometry 
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where 

D is the object radius and 
R corresponds to the Fourier frequency. 

The Radon Inversion Formulae [13 relates the data 
p(l,Q) to the function p(r,0) as follows? 

p(r,0) = ^ p(R,e) IRIdR de (7) 

Considering the radially synmetric flow into account 
p. no longer remains a function dependent on 0, so that we 
can choose any particular 0 and obtain the reconstructed 
p(r). For computational ease, and without loss of generality, 
we choose 0=0, Thus eqn,(7) beccsnes, 

li(r) = 7 / ptR) IRIdR d© (8) 

O -00 

where 

^R) = f p(l) dl (9) 

We see that p(l,©) has been replaced by p(l) since 
radial synmetry in p implies the independence of p with 
respect to the source position,©. 

Interchanging order of integration in eqn,(8), we 

get 

p (r) a °f ^R) ( 7 R^Cos© 

•“00 o 


m) |R 1 dR 
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Here integral over 9 is a Bessel Function [9]^ so 
that the inversion formula takes a particularly simple 
form, 

00 M 

H(r) = 211 / plR) Jo(2iirR)lRl c3R (10) 

o 

A practical implementation would require a finite 
cut-off for the variable, R, the Fourier frequency. Denoting 
this cut-off by R^, we get 

ii(r) 2* / p(R) J (2TTrR) IRl dR (11) 

“ o ° 

Keeping in mind the discrete nature of the data 
p(l), a reasonable choice of R^,, based on sampling theorem 
[10], is given by 

R^. = 1/2A1 

wrfiere, A 1 is spacing between the rays. It is to be noted 
that the expression given in eqn,(11) is similar to the 
point-spread function expression reported by Lewitt [s ]• 

In the design of reconstruction algorithm using 
Transform approach, it is often necessary to introduce a 
window function ^R) t 

W(R) = JW(R) IRl < 

io >RI £.Rc 
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giving eqn,( 11 ), the following form: 

li(r) ^ 2 % / p[R) J (2itrR) IRl W(R) dR (12) 

”■ o ° 

Two different types of window functions have been 
tested, viz., linear and sine functions, A sensitivity 
study has also been performed on the slope of the linear 
window function. The forms used are: 


Linear 




WCR) 



wile re 

0 < e ^ 1.0 


Sine 

Sin(itR/R^) 

W(R) = — 

(n R/Rc) 


2,3 Fan Beam Geometry 

In this case, one needs a rotating source and a 
ring of stationary detectors on the outer periphery. 
Scanning involves only motion of the source. 

In medical applications of CT, in order to reduce ' 
distortion due to movement of patient, data should be 
collected in shortest pxissible time. Ihis is facilitated by 
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using fan beam geometry of Radiation diverging from a 
source, passing through the patient and intercepted by 
detector elements. From the transmitted radiation, 
intensities are recorded while source moves in a circular 
path around the patient. 

The relationship between both the geometries [B], 
is as follows: 

1 = D Sin (T , © = p + cr 

p(l,e) = g((y,p) , 

2,4 T^yjjes, ,of^Aljg[qrithnis. 

For image reconstruction from projection data, the 
algorithms may be categorised as under; 

(a) Transform Methods [8 3 

(b) Series Expansion Method [tl] and 

(c) Distributional Approximation Method , 

(a) Transform Methods: 

Based on Radon Inversion Technique, the Transform 
methods basically rely upon analytic formulae . Transform 
methods [ 8 } ai"® broadly categorized as under; 

(i) Direct Fourier Inversion (DFI): in this 

meiiiod, direct Fourier Transform of projection 
data is taten. Subsequently 2-<i Fourier Inversion 
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leads to the reconstruction of unknown 
distributions, 

(ii) Convolution Back Projection (CBP): This 
method is based on convolving data with a 
suitable filter function. The back projection 
of this convolved data results in the recons- 
truction of unknown distribution. 


(j3) Se:^es Ejcpans^i^on Method* t 

In this method pixel-wise distribution of the 
function under-investigation is assumed. To obtain recons- 
truction of function in the region of interest [llj f 
suitable iterative and non-iterative methods are applied. 

The iterative SEf4s are: 

(i) Algebraic Reconstruction Technique and 
(ii) Simultaneous Iterative Reconstruction 
Technique etc. 

Whereas non-iterative SEMs are; 

(i) Angular Harmonic Decomposition and 
(ii) Polynomial Decomposition, 

( c ) DistrlJautjLonaJ^^ppr qx^i|[a^^^ 

The distributional approximation methods developed 
at IIT Kanpur involve assuming the distribution on the 
entire region of interest. 
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The RAP method [6,7 ] incorporates radial polynomials 
as the basis functions while the CSI Method [12] involves 
constant ring— values as the basis functions in its formu— 
lation.i 

The approach used in the present study involves the 
following steps; 

(1) Read data p in the form of a column vector 
corresponding to the spacing between rays, 

(2) Computation of p(R)* 

(3) Computation of p -values along the chord 

using eqns, (11,12) (involving Bessel Function), 

(4) Calibration of p -values to the density 
values , 

For the validation studies, data was generated (eqn,( 2) ), 
w*iereas in the eorperimental air-water flow studies, data is 
available from the detector counts. For computation of 
^R), a simple integration scheme was adopted (Appendix E), 

2,5 Two-rPhase Flow 

Sonte of the major problems in the nuclear industry 

are; 

(a) Measurement of irass-flow rate of two-phase 
flow, 

(b) Finding nature of f|x>w-around fuel bundle 
in Light Water Reactor (LMl), 



(c) Developing reliable two-phase flow measuring 
technique for various flow systems in a 
reactor. 

Earlier photographic methods for identification of 
flow regimes, using appropriate correlation to measure 
void-fraction were used. With the development of CAT, 
problem of determination of mass flow rate and identification 
of flow regimes has been simplified. 

Different two— phase flow regimes, namely, bubbly, 
slug, annular and dispersed flov\B can be investigated for 
determining void— fractions , However, to determine mass- 
flow rate, it is essential to know void-fraction without 
disturbing the flow field. This can be achieved by Tomographic 
Methods much more accurately than by other existing methods. 
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CHAPTER 3 
RESULTS 


• m. ^ h * 


3 .1 

We discuss the results for the simulated data in 
this section. We assume the object size/ pipe radius to 
be one unit and the ray spacing 0.1 , 

The algorithm was tested against the simulated 
data on the following symmetric distributions: 

li (r) = 1 .0 

p, (r) = r 

p ( r) - exp ( r) and 
p ( r) = exp (-r) 

The reconstruction is within + 2^ for all the 
test functions except in the peripheral region. The 
true and reconstructed values are shown in Fig, 2 , 

The algorithm has also been tested for annular 
flows, having distributions as under: 

0 


( r ) = 


1.0 1 

A A r 


= 0.0 


\i ( r) = 


0,0 
1.0 j 

0,0 j 


0 < r < 0,4 
0,4 < r < 1. 

0 <j r < 0,2 
0,2 < r < 0,4 
0.4 £ ^ < 1 .0 





a 



r - - - >0 

p (r) =1*0 



True value i 

l 

o Reconstructed! 

value ! 

{ 



r ► 

p (r)zexpk) 


I « 

t ? 

I t 

j ‘i 



Fig.2 Reconstructed results tor test-functions 



16 


In a caliberated sense, the above mentioned distri- 
butions represent the various types of density/void fraction 
distributions encountered in radially symmetric bubbly 
and annular flows. The reconstruction is showi in Figs ,3 
and 4, The solid lines represent the assumed functions and 
the dotted lines indicate the reconstructed values for the 
unfiltered and filtered versions; 

Testing of algorithm for Dirac Delta (S) function 
(at the center of the pipe) was also carried out. Both 
the filters viz, linear and sine were used for 6 —function. 
The results are encouraging, both for unfiltered and filtered 
versions (See Fig, 5), 

3.2 Results for Bubbly Air Water-Flows 

The test data is taken from the study of Ref, [3,4], 
Five different data sets for four different densities/void 
fractions were processed by this algorithm. The algorithm 
output, CTM, was calibrated to obtain the density. To do 
this, the data was processed for the following four known 
densities: 

Air, <P>- — 0,00 g/cm^, 

3 

Pine, < P>= 0,41 g/caea , 

Walnut;<^> =0,732 g/cts?, 

3 

Water, 1 ,00 g/cm , 

wtoere ” < > ” denotes the cross-sectionally averaged vala»s . 
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FIG. 3 RESULTS WITH FILTERED VERSION 
(ANNULAR FLOW). 
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FIG. 4 RESULTS WITH FILTERED VERSION 
(ANNULAR FLOW). 
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FIG. 5 RESULTS WITH FILTERED VERSION 
(DELTA FUNCTION). 



20 


The output CTN, was corrected to the base of air 
( CTN = 0,0) to eliminate the effect of the pipe— wall» made 
of plexiglass. We know, for air, the absorption coefficient 
(CTN) must be zero, whereas for air data set, we get non- 
zero values for the <CTNs> , This effect is due to atten- 
uation in the plexiglass pipe. In order to get the correct 
CTNs for all cases, we subtract the CTNs for the air 
(P= OpO g/cm ) case from the reconstructed values. This 
way, we can get correct CTNs and < CTNs > for any case. 

The CTNs obtained are listed in Appendix B, The calibration 
curve [<P> vs, ''iCTI'J> ] is shown in Fig, 6, by which, 
for any known <CTN>g the corresponding density value can 
be obtains do 

Similarly data for different scans were processed 
and vaJues of CTNs for the 12 rings in the pipe cross- 
section were achieved. For each radius, the p value 
obtained after calibration, have been plotted. Figures 7—11 
show the reconstructed density profiles for various 
densities, < P>, (or void-fractions <c^) for all the scans. 

In Fig, 12, the comparison of reconstructed densities 
with actual densities (measured by an alternate method [3,4]) 
have been shown. The various densities obtained by cali- 
bration nethod are marked to show the deviation of the 
reconstructed densities from the actual densities. The 
summary of <CTN> and<P>for all scare is given in 
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Appendix B, 

It is observed that the reconstruction is vdthin 
3 

+ 0,03 g/cm of the actual values. The variation of 
error vdth density has been plotted in Fig.13. The vari- 
ation of relative error with respect to density is shown 
in Fig, 14 and it ranges between + 3,57% to -6,67%, 

3,3 Results mth Filtered Version: 

The two filters viz, linear and sinc» applied to 
four known cases of distributions of 1,r,exp(r) and 
exp(-r), give encouraging results. Figures 15-18 show 
the distributions without filter and with these two 
filters separately, for these four known cases. The 
deviation of the reconstructed values from the true values 

is within + 0,039/ cm .( except for the perifSieral region) for 
all the filters tested in the study. It is to be noted that 
near the edges the density falls off sharply and a recons- 
truction incorporating a continuous function (Bessel function 
in this case) naturally leads to discrepancies in the peri- 
pheral region. However such errors do not significantly 
effect the reconstruction quality of the overall cross- 


sectional image. 
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LEGEND : 



FIG. 17 RESULTS WITH FILTERED VERSION. 
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CHAPTER 4 
DISCJ^SION 

The proposed algorithm, using the simplified Radon 
Inversion Formulae, has been demonstrated to be applicable 
to radially symmetric two-phase flows. The unfiltered 
version (EPS = 0,0), incorporating a unit step cut-off 
in the Fourier space, gives results comparable to the more 
general methods of the same type [13j, 

However, we observe that the unfiltered version is 
prone to slight overshoots/undershoots, vhile the filtered 
versions are relatively free of the aforesaid oscillationsi 
The results with the linear and sine filters show a reason- 
able Improvement in the reconstructed values (See Figs .15-18), 
These filtered versions are also comparable to the general 
tomographic methods involving the same filter-types [133, 

The bubbly flow projection data processed by this 
algorithm was in the range of 

. 3 

0,6 s: < P> 0,9 g/cm 

which is equivalent to, 

0*f ^ < a > < 0,4 , 

v^re g denotes the void fraction for the air water flow 
systems. The values of < Pqj> are within + 0,03 g/cm of 
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the true values, <P> This nagnitude of the deviation 
compares favourably compared to the other existing 
methods of measuring density/void fraction in two phase 
flows . 

Vi?e also note that the results reported here are 
comparable to the other algorithms developed [7,12] for 
radially symmetric situations. 
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APPENDIX B 
TABLE -B1 

(Results for Air, Pine, Walnut and Water) 


3 

Case < p > » g/cm <CTN> 

Air 0*000 0,000 

Pine 0,410 0,081 

Walnut 0,732 0,146 

Water 1,000 0,212 






TABLE-B2 

(Reconstructed Densities) 





3 

4 

5 


0,57 

0,56 

0.59 


0.7 


O 48 


0.9 


r '•fc . \ a 


0,^72 

0,70 

0.67 

0,69 

0;72 


0.82 

0,80 

0,78 

0,81 

O482 


0,91 

0.90 

0,89 

0.92 

0i91 




,:f4 ■'■WtC- ■- % 


I m ^ 


TABLE -B3 


(Reconstructed <CTNs>) 


^\<P >g/cm^ 

'-•n n _■%. ■^-. ^ _ 

V A- \ ^ 



Scan 

’^oT"' X 

0,6 

0.7 

0*8 

♦ •v .»A K_ ^4 4 ,J, -%. 

0i9 

1 

0,1165 

0.1443 

0,1638 

0,1829 

2 

0,1150 

0,1383 

0,1601 

0,1808 

3 

0,1114 

0.1310 

0,1572 

0,1792 

4 

0,1097 

0,;1370 

Oil 610 

0.1847 

5 

0.1157 

0J445 

0,1653 

0.1843 






■%. : - %-■% 



1 


RESUJuTS for fun = 
p(L) = DATA 


WITHOUT 

FILTERED 

VERSION 



li— A T’ 

PCT.) 

R'^.MnLT. 

FCR) 

TNT-VAT, FRPQP 

O.^OO 

0.200 

0.300 

0.40^ 

0*50O 

O. 60 O 

n iQfi 

o'PoO 

O.OQf’ 

r.oo^ 

2 . non 

j OQfi 

1 Ins'-* 

1 , 0 -jP 

1 .R 33 
1.737 

1 . Aon 

1 .4?5 

1 ,900 

n. P 7 ? 

0 ^ n 6 r? 

n.ooo 
n. 1 on 

0 0 AO 

n " 7 0 n 
n>on 
n.Ajn 
f'.AuO 
n.7jn 
n , P u n 
n ^ 0 n 

1 loun 

1 .f'oo 

l.non 

1.^00 

1.030 

1 .non 

1 .non 
l.non 

1.000 

1 .non 

1 .non 

1 .000 

- 0:033052 

0.00455A 
i*2i7^4n -0.01734Q 
0.004357 
1.01773? -0.017737 
1.305409 -O.OOS^oQ 
1*02^52^ -0-005603 
1.040427 -!), 040427 
’->-011^75 
1 . 16741 ? -0.167417 
0.27041P 0 . 72953 ? 


WITH 

LINEAR FILTER 




L-.VAL 

PC L) 


TCP) 

I.vr-VAL 

E8R3R 

3.>.'30 

0,1 30 
3. -230 
0.3 30 
3,430 
3.-530 
3,630 
0.-7 30 
3.830 

0. 930 

1. -.OO 

2, U3D 

1 .'PSO 

1.732 

l:| 2 i 

klU 

3 , -3 3 3 

3. -3 DO 

0 .-no 
0.-23 0 

i-.m 

1-Mi 

0.-7 33 

0,-3 00 

0. -9 3 0 

1. -030 

1.'333 

1 .'33 3 
1.'30 3 
.033 
1.033 

1 .033 
1.033 

1 ,030 

1 .03 3 

1 .-333 

1 .033 

0.982484 
0.978126 
0,979857 
0.978213 
0.9771 22 
^977493 

0,372651 

0.017515 
0.021874 
0,0231 43 
0.321793 
0.322978 

i-Mmi 

3.025243 

0.027874 


WITH 

SING FILTER 




l-val 

PC 00 

FCr) 

i«r’“VAL- 





2 * 

1.393 
1 . 9 5-3 
1.333 
1.333 
1.7 32 
1.&33 
l.i?'3 
1.233 
0.372 
0.333 


03 3 
-i 03 
■203 
033 
•133 
3. '5 3 3 
3.|03 
3.'/33 
3.033 


3. 

3. 

3. 

3 . 


3. '3 3 3 
1.033 


I'.OOO 
l.ODO 
1 .033 
1.033 
1 .0 30 
1.033 
1.033 
1.033 
1.033 
1 .'333 
1 .0,33 


1.309075 

I. 034987 

J. 037303 

1.0369J4 

1,037557 
1 .010953 
3 .010397 
1:019755 

1.043919 

0,996265 

0.369192 


- 0,039375 

- 0.034987 
-O.O 373 O 8 
- 0,03691 4 

- 0.037557 

- 0.013953 
“ 2.^12397 
- 0.019753 
- 0.04391 9 
0.003734 
0.631838 



RESULTS FOR FUN 
P(L) *= DATA 


WITHOUT 

l-val 

1.000 
. .100 
; ,200 
K.300 
i.4O0 
^ ,i>00 
,0.600 
4.. 700 , 

. .800 

. .90o 

1 .000 


FILTERED 

PCD 

1 .OOo 
1.025 
1.07i 
1.123 
1.167 

1.195 
1.1°6 
1.153 
1.044 
0.814 
O.OOj 


VERSION 

Rs*Hr,L 

O.OOJ 
O.inO 
0.200 - 
O.30j 
0.40i- 
0.50U 
0.60 0 
0,700 
0 . 80 ,? 
0.90L, 

1 . 'J 0 


F(H) 

0.000 
0.100 
0.2 00 
0.300 
0.400 
0,500 
O.oOj 
0.700 
o.soo 

0.9 00' 
l.OO'O 


INr-v'AL 

0.078407 
0.l03b94 
0.221228 
0.298952 
0, 419104 
0.507586 
0.60655U 
0.741505 
0,789629 
1.066993 
0,254b35 


ERROh 

-0.078407 
-0.003o94 
- 0.021228 
0. 001046 
-0,019104 
-0.007586 
-0.006550 
-0.041505 
0.0i037i 
-0.166993 
0.745365 


WITH LINEAR FILTER 


6-VAL 

: .000 
‘ .100 
;.20o 

'>.300 
; .400 
*'.500 
: .600 
V.700 
^.800 
; .900 

i.ooo 


P( 60 

1.07c 

1.025 

1.071 

1. 12i 

1.157 

1.195 

1.196 
1.153 
1.044 
0. 3H 
D.OOO 


O.'OOj 
0.400 
0.200 
o.'3o: 
0.0 oo 
0.500 
0 • 5 0 0 
D.OOO 
0.'3 0D 
0.900 

1.000 


F(R) 

o.too 

0.100 

0.20D 

O.SOO 

*^.90 V 

0.500 
0 . 6 0 J 

0,700 
0.600 
0,900 
i .UOO 


I^Ir-VAL 

0.12^132 

0,157358 

0.236552 

0.322135 

0.411923 
0.5D4S?? 
5:*553349 
0, 687013 
0.779142 
0.803276 
0.338295 


iFROR 

-0.'l24i32 

-0.057350 

-0.036552 

-0.022136 

-O.0U92O 

-0.004330 

0.035651 

0.020353 

0.096724 

0.661705 


WITH SINC FILTER 


L-VAL 

0.000 
*’.100 
; .200 
.*.300 
: .400 
f^.500 
D'.&OO 
e.70o 
7.800 
: ,900 
1.000 


P( LO 

1.000 
1.0 25 

1.123 

1,157 

1.195 

1.196 
1.15 3 
1.044 
0.314 
O.OOO 


R'SHAL'Li 

0.000 
O.'iOo 
D,'2 0 0 
0.000 
O.'JOO 
3.500 

l-Ml 

O.'BOO 

0. 900 

1. 


FCr) 

0,000 
0.1 00 
0.200 
O.SOO 
0.400 
0.500 
0.500 
0.700 

o.eoo 

0,900 

1.000 


iWr-VAL 

0.036652 
<^.1 25537 
M1I543 

0.4i?496 

0.514872 

0.614153 

0.722440 

0,846333 

0.896221 

0.343932 


iPROR 

-0.035552 

-0,025537 

-0.013544 
-0.012495 
-0.014372 
-0.014153 
-0.022440 
-0,0463 30 
0.003779 
0.656053 



results 

FOR FUN 

= EOCP(it) 



P(L) = DATA 




WITnoUT 

F ILT£.kii.L 

V£,RS10W 



b-y^AL. 

PCD 

RSMaLu 

FCft) 

INI-VAL ERnOn 

,^.000 

'.lOQ 

.200 

:^.3O0 

...400 

>. 50 v; 

, .600 
./. 7 O 0 

.U800 

»*.yoo 

1,000 

3.437 

3.457 

3.491 

3.513 

3.5’7 

3.406 

3.494 

3 . 21 s 

2. 869 
2.219 

0.0 0., 

u* . 0 0 

0 . i 0 
0.20.- 
D. 3^>V 

0. iOO 

0.50 0 

0 . t» 9 0 
0.700 

O.rfOj 

0.900 

1.000 

l.ODO 

1.105 

1.221 

1.350 

1.492 

1.649 

1.822 

2.014 

2.22b 

2.460 

2.718 

1.138516 -0.138516 
l.l 0206 j 0.0D31I1 

1.274623 -0.053220 
1,343425 0,006434 

1,543105 -0.051200 
1.668565 -0.019844 
1.840291 -0.018172 
2.127485 -0.113733 
2.198582 0.02637^^ 

“2.455770 
0.693351 2.024331 


WITH 

L-VAL 


vAU 

f.3O0 
:• .400 
'•500 
'.600 
: .700 
. ,800 
: .900 
1,000 


LINEAR FILTER. 


PC LO 

R'SMAL'U' 

FCR) 

iNr-VAu 

3.437 

0.000 

1.000 

1.130545 

3,457 

0.10 0 

1 .1 05 

1.151727 

3.491 

0.200 

i.2h 

1,250079 

3.513 

0.300 

1.050 

1.371043 

3.527 

0.'4 0 0 

l.‘492 

1.500060 

3.495 

0.5 0 0 

1.649 

I. 648385 

3.43 4 

0.5 0 0 

1 .S22 

1.739902 

3.215 

0.'7 0 0 

2.014 

1.933251 

2.8S3 

0.500 

2 ,22b 

2.176438 

2.213 

0.9 9 0 

2.450 

2,20631 9 

0.09 :• 

1.000 

2.718 

0.927383 


-2«U?545 
• 0, 036565 
•0,038576 
- 0.021104 

- 0.000164 

O.O 222 I 7 

0.030491 

0.049103 

0.262694 

1.730599 


WITH SING FILTER 


li-tfAL 

P ( uO 

RSWh 

FCR) 

>.D00 
0.100 
f’.2O0 
» , 3 0 0 
0!.40D 
?,500 
1>.600 
: .700 
0,800 
'.900 

1 ,000 

3.437 

3.437 
3,491 
3.519 

3 - 5 77 
3.495 

3. -4 0 4 

3.215 

2.359 

2.219 

0,000 

0.000 

0.100 

0.'200 

0.>3O0 

i-.ni 

1.00 0 

1.000 

1 .'350 
1.'492 
1.649 
1.622 
2.014 
2.226 
2.160 
2.718 


1%’r^VAL 

1.106184 

1.144564 

!;256031 

11330293 

1.522628 

r.53Pl38 

i;85l260 

2.076739 

2.354704 
2 453574 

0’93748& 


£ 8 R 0 S 

-0.'l06l9i 

— 0 ,0 34 528 
-0. 039374 
- 0.039304 
- 0.039417 

-0.039i4l 

-0.052395 

-0,-l29l54 

0.006023 

!.. 730795 



RESULTS FOR FUN = exp(-r) 
P(L) = DATA 


WITHOUT 

filtered 

VERSION 



Ij-VAl 

UCb) 


F(k) 

INf-tfUb ERROk 

oou 
.100 
=. ,200 

1, • J 0 u 

'. 4 OU 

...500 

.60u 

. .700 
.800 
;.90 3 
l.OOO 

1.264 

1.234 

1.159 

1 . 0? 4 
0.98o 

0.87g 

0.754 

0 . 0 4 
0.504 
0.3?3 
O.OOu 

0 . U 0 'J 

3 . 1 0 J 

0.203 

0 . 301 ; 

0. 40v> 

0 • 5 0 G 
3.60 3 
0,70o 
o.Bo: ■ 
0.900 
1,003 

1 .OOQ 
0.905 
0.819 
0.741 
0.676 
O. 6 O 7 

0.549 

C .497 

0,449 

0.407 

0,368 

S-?22768 “0.004337 
n'lUiU S*2^358i 

0,676513 -9,006i9i 
0.607694 -O.OOlisI 
$•1^^756 - 0.031944 
£«|lU|^ -0.014774 

0,444774 0,034555 

0. 468666 -0,062096 
0,105423 0.26245b 


WITH 

LINEAR FILTER 




L-VAL 

PC 50 

8sMA,LLi 

FCR) 

191-VAL 

SRR3R 

'...DOO 

1^.100 

%200 

V' . 300 
. .400 
.500 

V .600 

» .700 

* • 8 0 0 
.900 
1.000 

1.254 

1.234 

1.159 

1.33? 

0.996 

0*379 

0.751 

0.534 

0.343 

D . g 3 

D . '1 0 3 
j.'2 0 3 

3. -30 3 

3,4 0 3 

3. '5 03 

3. '603 

3. '7 03 

3. '3 0 3 
3.903 

1 .903 

1,000 

3.935 

3.919 

3.741 

3.673 

3.637 

3.549 

3.497 

3.449 

3.437 

3.'3 5o 

0.874447 

§: 

0.536635 

O.53U10 

S-f32955 
0. 43^275 
0.375638 
0.151383 

olpilp 

0.043012 

5.032259 

0.019396 

0.017501 

0.011954 

0,030931 

0.216799 


WITH 

SINC 

FIjuTj:,R 



o-val 

PC 

u7 

RS'IaL'O' 

FCr) 

. .300 

1. 


3,033 

1.030 

. . 1 0 0 

1. 

234 

3. -1.3 3 

0.935 

‘.200 

1 • 

159 

3. '2 3 3 

3.619 

* * 3 0 0 

1 . 

384 

3.030 

3.741 

G . 400 

0. 

985 

3. *4 0 0 

0.670 

. 5 0 0 

0. 


5. '5 0 0 

3.637 

■ . 60 0 

0 • 

754 

3. '5 0 0 

3.649 


0. 

■§|o 

3 . '7 3 0 

3, <497 

• B 0 5 

0. 

0. 

534 ■ 
343 

3 '! 3 3 

3.449 

3.437 

t • 0 0 0 

0. 

330 

1.000 

3.358 


iNr-VAL 

0.929526 
0,fi8S792 
o:«i 4566 
0.739895 
0.671 394 
0.639357 
Ol552ll5 
0,533478 
0.455285 
0.435823 
0.1 44746 


£RR3R 

0.070474 

0.004155 

0,000923 

-0.030774 

-0.032826 

-0.033304 

-0,035393 

-0.015955 

0.030647 

0.223133 



1 - iJ, 

u n T . ? ,% ■'* j •” ii r i i i r v j i , t r !,■ 


A!( 

1 > A A 


F 

f K V 

! 1 ” 

1 r 1 1 

0 . 0 0 y 

1 

,uOt) 

1 ' ■ ' 

' , /9s 

0.1 

1 

.CA(? 

f t 

' .s'^l 

i ■ . 2 ^ () 

1 

,(>aC 

> 

i .SOU 

U . 3tSt 

1 

• M 0‘.i 



0 A A u 

Cl 

, 0 A u 


* ! 4 

• ' ' V 

oI^om 

1 ; 

.900 


1 # 

c.cro 

0 

• VOO 

y ( 

1 ,f r,,i 

r.70y 

fi 

. or . i ) 

F "" 


C.i:ro 

f'i 

‘.C.OC 

4 \ f > « . 

( , (■'<( 

(■•*90(' 

0 

'.ono 


u, i, ''C 

1 .COO 

(1 

',yoo 


INT-V'M, 

0.^7^2'?2 

i.og.'ist?’/ 

0, q7'lS04 

1 , UiPg'vS 

-0,0215!^'} 
-O.OQ'iTOo 
*■> 0i?57(j 

-0,0l2s}'’y 

-'■>.00 147 ‘I 


error 
0 . 12S708 

- 0,094507 

0,025306 

'0,189855 

764270 
0.02! 5 R 8 
0.0(1 57 Of, 
-o.Oj 7578 
0,01 7079 

-O.fKiR^nv 

0.003404 


'■•cr.in, ’’P. 


fiTtlFAR FltiTEP 


Cl A 1 / 

D A i’ A 


t', iitS- 


( . /''s 

? C'M 1 

' . f''* ! 


i:<.s7u 

4 ^ ^ 


*3 0 1 . 

( , ' fu 


(.,r r , 

1 

(’.( A' 

F' ♦ * t " 

t , i'('' 

4 ( ' i 

' . 1 . o« 


A . C 0 (. 


RS'*AT.L FfR) 

O.OOO 1,(if>0 

O.tOO 1,000 

0.200 l.OOti 

0,300 3.000; 

O .400 0,000 

0,500 0,0 00 

O.f.oo 0,0' AO 

0.700 O.OOO 

0.000 

0,000 0 ,on(, 

1 .COO O.V'O 


IHT-yftu error 

0,950976 0.049074 

0,960556 0.039444 

0.959434 0,041555 

0,886079 O.moSi 
0.140691 -0.340691 

0,034905 - 0,034996 

0.0394f;7 -A, 019457 

0,00795" -9,007950 
9,005277 -A, 005277 
0,903779 -O.no3778 
0.007179 -0. 007179 


«-l.r.lJl/J'S *''jTh 51 or Flf.TFH 


1,-V7l. 


(..r.Of 

(., ir v> 

C, If <- 

U * of'c 
( ,lf'0 
i 

( .f'f'O 
3 .t Cf 


b A T A 
C . bf^C 

(',775 

(■■.503 

.57*1 
.1 <'( 
.( tH 

.( f'A 
.(<('<• 

A 0(1 

.C0(, 


RKMAT.L 

0 . OOn 
(.',100 
(> , 200 
0,300 
0.400 
0.500 
0,600 
0.700 
O.BOt 

0,900 
1 .000 


FftO 

1 ,00(' 
3 .COO 
1 .(.00 
1,000 
0.000 
0,0 00 
O.OOO 
0,000 
0,09 0 
O.OOO 

c.coc 


XNT-VAb 

1 .021504 
f .037216 
3,056784 

0,990315 

0 . 15090 b 

.0.007565 

0.007971 

.0,000975 

0.000090 

0.000230 

.0,000278 


ERROR 

-0,021504 
-0.037216 
-0,056784 
0,001665 
-0,150906 
0.097565 
-A, 002971 
0,000975 
-0,000096 
-0,000230 
0,000278 




P-UMctxpk't 0 , 1. 0, 



pe^ijT. 

TR WITHOUT 

ETLTER 

r,-v4i., 

0 A r " 

RS'^ATl. 

FfPl 

0,H9() 

U, 40u 

olooo 

O.UOO 

0 , 1 0 0 

0 . 42 r 

0,1 Co 

0 , 0 0 0 

0,2<^0 

0.69 3 

0,29a 

1 .000 

0, 

0.529 

0,300 

1.000 

y , 4 0 y 

0 , 0 9f) 

0.400 

0.000 

0,590 

0 , 0 9 n 

0.500 

0,090 

y.e^o 

0 , 9 9 0 

0,6C0 

0,000 

f'.TOn 

0 . 0 0 i) 

O. 7 OO 

0.000 

0,B9,) 

0.0 9.J 

O.ROo 

0,000 

0,990 

0, yoo 

0,900 

0.000 

1 .09'» 

0,«9() 

i.eoo 

0,000 


IMT-VAL 


0 . 004273 
“0.1 682'>0 


0,757535 


1.180902 

0.277719 

“ 0.032493 

“ 0,000701 

0,012591 

“ 0,015140 

0.012105 


- 0.005490 


error 

“ 0,004273 

0 .U 0250 

0,247465 

“0,189902 

“ 0,277719 

0,032493 

0.000701 

“ 0,012591 

0,015140 

- 0.012195 

0,006490 


’’eSlJT.TR WITH linear K'lr.TER 


b*V A fj 

D A r A 

RS'*ALl. 

FfRl 

INT-VAL 

error 

0 • f.) 9 0 

0,40,, 

o'. 000 

0,000 

-0,031784 

0.031784 

0,10 y 

0.42fj 

0.100 

0,000 

0.1 11217 

-0,111237 

0,29u 

0.691 

0.290 

1,000 

0,662904 

0. 117096 

0.390 

0.529 

0.300 

1,000 

0,865200 

0,134800 

0 . 400 

C'.OOO 

0,490 

0.000 

0,329495 

-0,329495 

0 , RPd 

O.OOo 

0.590 

0,000 

0.030962 

-0.030962 

O.hOo 

o.y9f> 

0,690 

0,000 

9.017311 

-0.017311 

y,70ti 

y • oOo 

0.700 

0.000 

0.006223 

-0,006223 

0. H90 

o.yoo 

0.800 

0.000 

0,004505 

-O.O04505 

0, wOo 

0,0 Oo 

0,900 

0,000 

0.00299? 

-9,002997 

t . yoi) 

C'.yoo 

Ilooo 

0,000 

0.001633 

“0,001633 


ReSUr.TS WITH STNC PILTBR 


L-VR L 

narR 

RSWRLL 

fcri 

INT-VRL 

ERROR 

y.yoo 

0. 100 
0,200 
0.300 
0,400 

0,500 

0.600 

0,700 

0 ,H 90 

0,900 

1 ,000 

0,400 

€.428 

0.529 
( . ooo 

O.OOo 

C .090 

0,990 

O.nOn 

0.000 

C . OOo 

O 10 OO 

0.100 

0,200 

0,300 

0,400 

0,500 

0,600 

0,700 

0,800 

0,900 

1.000 

0,000 

0.000 

1,000 

1,000 

0,000 

0,000 

oIo|S 

ttn 

0.000 
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